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''I.  INTRODUCTION 


"With  the  increasing  importance  of  lasers  In  our  technological  world 
and  the  continually  increasing  power  levels  becoming  available,  it  is 


inevitable  that  laser  weapons  will  play  an  important  role  in  various 


phases  of  modern  warfare. x Any  alert  observer  of  the  news  media  understands 


that  lasers  have  been  and  dpntinue  to  be  developed  for  a variety  of  military 


purposes  among  which  are  offensive  weapons,  defensive  weapons  (antiaircraft. 


antipersonnel,  etc.),  guidance  of  bombs  and  missiles,  and  range  finding  and 


other  types  of  reconnaissance. -> As  dependence  on  lasers  increases,  there  is 


a need  to  know  how  these  intense  light  beams  are  affected  by  airborne 


particles,  either  intentionally  or  accidentally  produced,  whether  man-made 


or  natural.  Obvious  defensive  mechanisms  against  laser. weapons  might  include 


smoke  screens,  i.e.,  high  concentrations  of  particles  selected  to  effectively 


stop  certain  laser  wavelengths.  Against  this  background  of  laser  development. 


this  study  has  begun  to  examine  the  processes  by  which  laser  light  is  atten- 


uated through  interaction  with  a collection  of  small  particles. 


A.  Purpose  of  This  Work 


The  direction  of  this  work  has  been  toward  determining  the  most  effective 


type  of  particles  for  the  purpose  of  attenuating  light  from  CO2  lasers.  These 
are  the  most  common  types  of  high  powered  lasers  radiating  near  9.6  and  10.6  pm 


in  the  infrared.  Since  there  are  so  many  kinds  of  particles  that  might  be  used, 


this  wor k has  focused  on  only  two  kinds  of  particulate  systems.  They  are: 


(1)  silicate  smokes  — very  small  ( .05  pm)  silicate  particles  and 

(2)  glass  bubbles  — hollow,  thin-walled  glass  spheres  ( 10  pm  o.d). 
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These  were  considered  to  be  good  candidates  for  producing  appreciable 
attenuation  of  the  CO2  laser  wavelengths  for  two  different  reasons 
relating  to  the  processes  involved  in  stopping  a beam- of  light. 

These  two  processes  are  scattering  of  light  out  of  the  beam  and 
absorption  of  light  within  the  particles,  the  sum  of  these  stopping 
processes  being  known  as  extinction. 

extinction  = absorption  + scattering 

(Mathematical  expressions  for  these  quantities  will  be  given  in 
Section  II.)  The  glass  bubbles,  already  available  commercially,  were 
expected  to  give  large  extinction  by  the  scattering  process  since  the 
effective  sizes  are  approximately  equal  to  the  infrared  laser  wavelengths. 
Because  of  the  thin-walled  nature  of  the  bubbles,  the  mass  per  unit  of 
extinction  would  be  kept  low;  that  is,  the  bubbles  were  expected  to 
provide  high  extinction  per  unit  mass  by  way  of  their  scattering 
properties.  The  silicate  smoke  particles,  on  the  other  hand,  were  known 
to  have  quite  high  absorption  efficiencies  in  the  vicinity  of  10  pm 
wavelength. 

For  the  purpose  of  determining  the  infrared  extinction  from  the  two 
types  of  particles,  two  approaches  may  be  used.  The  extinction  may  be 
calculated  from  known  optical  properties  of  the  solids,  or  the  atten- 
uation of  light  through  a collection  of  the  particles  can  be  measured 
directly.  One  advantage  of  the  calculations  is  that  other  parameters 
such  as  size  distribution  can  be  easily  changed  in  the  calculations, 
while  the  experimental  technique  is  limited  to  the  specific  sizes  of 
particles  available  or  produceable.  The  calculations,  however,  are 
fairly  difficult  and  input  parameters  such  as  optical  properties  and 
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appropriate  size  distributions  must  be  experimentally  determined. 

Also,  the  calculations  can  only  be  accomplished  for  ideal  geometries 
such  as  spheres  or  spherical  shells.  The  applicability  of  these 
calculations  to  other  shapes  needs  to  be  examined  by  comparison  with 
measurements. 

The  plan  of  this  work  has  been  to  determine  optical  constants 
and  size  distributions  necessary  for  doing  theoretical  calculations  of 
extinction  for  both  kinds  of  particles.  Extinction  measurements  are 
made  to  determine  the  validity  of  the  calculations.  Recommendations 
can  then  be  made  based  on  the  comparison  of  results  for  the  two  types 
of  particles.  Finally,  the  most  promising  parameters  can  then  be 
predicted  from  the  calculations. 

B.  Limitations  of  this  Work 

This  work  will  be  restricted  to  the  two  representative  types  of 
particles  — glass  bubbles  and  one  of  the  most  promising  of  the  silicate 
smokes,  olivine  smoke.  Further  work  on  solids  of  different  chemical 
composition  will  be  left  for  other  investigations.  A further  limitation 
is  that  we  do  not  study  the  problem  of  how  the  particles  are  to  be 
produced  and  dispersed  as  an  aerosol.  These  may  be  important  and 
difficult  problems  of  a practical  nature,  especially  since  the  extinction 
is  known  to  be  somewhat  dependent  on  the  degree  of  agglomeration  of  the 
dispersed  particles. 

C.  Plan  of  the  Presentation 

This  report  will  summarize  the  research  methods  used  and  the  results 
obtained  during  the  past  year  under  the  support  of  NASC  contract 
N00019-77-C-0293. 


The  next  section  (II)  will  surnmarize  the  necessary  theory  by  way 
of  the  main  equations  used  in  the  scattering  theory  for  spherical 
particles  and  spherical  shells,  and  theory  needed  to  analyze  specular 
reflectance  data  from  solids  to  determine  complex  optical  constants. 

An  exhaustive  treatment  of  this  theory  is  not  necessary  here  since  it 
is  covered  in  various  readily-avai Table  references.  The  summary  of  the 
main  equations  is  desirable,  however,  in  order  to  clearly  define  the 
terms  used  and  to  explain  the  methods.  Following  the  theory  section, 
results  of  measurements  and  calcualtions  on  each  of  the  two  particulate 
systems,  glass  spheres  and  silicate  smokes,  will  be  presented  in  turn 
(Sections  III  and  IV  respectively).  Each  of  these  sections  will 
include  the  necessary  discussions  of  optical  constants  determinations, 
particle  size  distributions,  extinction  calculations,  and  extinction 
measurements.  Section  V will  then  compare  the  resulting  extinction 
spectra  for  the  two  particulate  systems,  and  will  present  conclu- 
sions based  on  this  work. 


II.  THEORY 


A.  Theory  of  Extinction 

A well  defined  beam  of  light,  such  as  a laser  beam,  is  attenuated 
as  it  passes  through  a collection  of  particles.  The  attenuation  is 
caused  by  two  processes  — scattering  away  from  the  original  direction 
by  the  particles,  and  absorption  of  light  within  the  particles  with 
subsequent  conversion  of  the  photon  energy  into  heat.  The  transmission 
I/I0  of  a beam  of  light  having  initial  intensity  I and  final  intensity 
I after  traversing  a distance  & through  a collection  of  uniformly  dispersed, 
randomly  oriented  particles  can  be  written  as  follows: 


I/I0  = exp  (-ftoext1  * • 


(il-1) 


HA*. 
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N Is  the  number  of  particles  per  unit  volume  in  the  light  path  and  a 


ext 


is  the  cross  section  for  the  extinction  process,  which  can  be  separated 
into  the  separate  scattering  and  absorption  parts 


ext 


sea 


abs 


The  major  computational  problem  in  calculating  extinction  by  particles 
is  to  determine  the  extinction  cross  sections,  oext  . In  general  this 
involves  a solution  to  the  electromagnetic  boundary  value  problem,  which 
is  quite  untenable  in  the  cases  of  particles  having  arbitrary  shapes  and 
arbitrary  sizes.  The  solution  to  this  problem  for  the  special  case  of 
spherical  particles  was  developed  early  in  the  present  century  by  Mie  (1908) 
and  by  Debye  (1909),  generally  referred  to  as  Mie  theory.  The  equations  for 
their  solution  are  complicated  but  are  nicely  given  in  a number  of  more 
modern  references  (Stratton  1941,  van  de  Hulst  1957,  Kerker  1969).  The 
necessary  input  parameters  are 

the  wavelength  of  the  light  - X , 
the  particle  radius  - r, 
and  the  complex  index  of  refraction  for  the 
solid  at  the  wavelength  - m = n + ik. 


In  recent  years  high  speed  digital  computer  techniques  have  permitted  the 
Mie  equations  to  be  readily  solved.  In  this  work  a Fortran  program  based 
on  work  by  Dave  (1968)  has  been  used. 

The  boundary  value  problem  for  concentric  spheres  was  solved  much  later 
(Guttler  1952).  The  equations  are  more  complicated  and  consequently  the 
computer  program  is  more  difficult,  but  similar  solutions  exist.  The  required 
input  parameters  are  similar  except  of  course  that  both  inner  and  outer  sphere 
radii  must  be  included,  and  separate  optical  constants  must  be  supplied  for 
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the  inner  sphere  and  the  outer  sphere.  In  this  work,  the  inner  sphere  was 
always  considered  to  be  empty  space  ( m -•  1.0  + iO.O),  and  the  outer  shell 
v/as  glass,  with  optical  constants  measured  in  this  v/ork. 

For  non-spherical  particles  there  are  no  general  solutions  for  arbitrary- 
sized  particles.  A sometimes  useful  approximation  can  be  obtained  for 
ellipsoidal  particles  in  the  so-called  Rayleigh  limit  (where  the  size  of 
the  particle  and  its  optical  dimensions,  n x length,  are  very  small  compared 
to  the  wavelength).  Using  this  approximation  one  can  sometimes  get  insights 
into  the  effect  of  varying  shapes.  Usually,  one  assumes  that  the  spherical 
approximation  applies  reasonably  well  to  non-spherical  particles  also,  but  it 
is  well  to  continually  remember  that  this  assumption  is  being  made.  Reference 
to  experimental  measurements  on  real  (frequently  non-spherical)  particles  is 
then  a good  idea  for  delineating  the  limits  of  validity,  of  the  simplifying 
approximation.  This  is  the  approach  we  have  used  in  this  work.  All  calcu- 
lations presented  are  for  the  spherically  symmetric  cases  — spheres  and 
concentric  spheres. 

For  the  cases  of  spheres  or  spherical  shells,  the  light  beam  transmission 
through  a path  length  l can  be  written. 

I/IQ  = exp  ( -NQext  tv  r2)  & 


= exp  (-  £ QOYt  tt  r2)£ 


(II -2) 

- J-\0. 

V wext 

where  Q . = Q„  „ + Q_._  . The  Q's  are  known  as  efficiency  factors  for 

GXt  SCa  aDS 

spheres  which  are  obviously  transformed  into  cross  sections  by  multiplication 
by  the  cross  sectional  area  of  a sphere,  irr  . N is  the  number  of  particles 
per  unit  volume,  n/V.  Although  this  expression  is  quite  acceptable  -for 
calculations,  it  is  sometimes  useful  to  relate  the  transmission  to  the  mass 
concentration  of  particles  c (mass  of  particulates  per  unit  volume  of  space). 
For  a monodisperse  (single  size)  collection  of  particles. 


- ir-~  ■ 

~ " -• 


m _ n 4/3  u r p 

c ‘ T ~ ~ 


(H-3) 


where 


n = number  of  particles  in  volume  V 
r = particle  radius 

p = mass  density  of  particulate  material 


It  follows  simply  that 


JL  - 
V 


4/3  7T  r3  p 


which  is  put  into  Eq.  ( I 1-2)  to  yield 


1/1  = exp  ( 2 Qn  /)i  (II -4) 

0 4/3  Tir3 


= exp  ( - ap)  ~£ 
= exp  ( - op)  £ 


( 1 1-5 ) 
( I 1-6 ) 


in  which  a is  the  mass/area  of  particles. 

The  extinction  cross  section  per  unit  volume  of  the  solid,  aD,  is 
defined  by  Eq.  ( I 1-5)  and  is  related  to  the  extinction  efficiency  factor 
for  spheres,  Qext,  by  comparison  with  ( 1 1-4)  as  follows: 

a„  = Q irr2  , = 3/4  ^ ext  . (1 1-7) 

p 4/3  n r°  r 


Many  authors  prefer  to  calculate  extinction  cross  sections  normalized  per 
unit  mass  of  the  solid  material,  a,  sometimes  referred  as  the  mass  extinction 
coefficient.  This  is  simply  related  to  our  ap  by  the  solid  density, 

“ - = 3/4  -L  ■ (II-8) 

with  the  corresponding  transmission  equation  being 

I/I  = exp  (-  ac£  ).  ( II-9) 


I'r'i'.f.'f1'  i~iv~ 


By  measuring  the  mass  density  of  particles  collected  on  a filter  or 
substrate  (a  ) and  the  transmission  through  the  coated  substrate  (I/IQ), 
and  using  Eq.  ( I 1-6) , the  ap's  can  be  determined  experimentally,  leading  to 
an  experimental  determination  of  Q , for  spheres  from  Eq.  (I 1-6).  These 

CA  t* 

equations  provide  the  connection  between  theory  and  experiment  for  mono- 
dispersed  collections  of  particles. 

If  there  is  a size  distribution  of  spheres,  the  volume-normalized 
extinction  cross  section  is  determined  from  the  particle  size  distribution 
function  N (r)  as 


r 2 
a p = N (r)Q(r)TT  r dr 

N (r)  4/3tt  r3  dr 


(11-10) 


rather  than  Eq.  (I 1-7).  N (r)  is  the  number  per  unit  volume  of  particles 
having  radii  between  r and  r + dr. 

In  the  case  of  speherical  shells,  once  the  Q's  have  been  calculated,  the 
Op's  are  determined  through  division  by  the  volume  of  the  shell  (rather  than 
of  the  sphere)  or  by  the  volume  distribution  of  the  shells  as  follows: 


4/3  *(r3  _ r.3) 
o i • 


4/3irr03  [1  - (r^)3] 


(11-11) 


[N  (r0)Q(r0)7rr0Z  drQ 
J « (rQ)4/37r rQ3  [1-  (rQ/r^  )3]drQ  , 


(11-12) 


where  rQ  and  r-  are  the  outer  and  inner  radii  of  the  shell  respectively.  Note 
that  the  radius  distribution  in  Eq.  (II  12)  is  over  outer  radii.  In  this  work 
we  consider  the  ratio  of  inner  to  outer  radius  to  be  constant. 
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B.  Theory  of  Optical  Constants  Determinations 

The  non-coherent,  linear  optical  properties  of  a solid  are  summarized 
by  the  index  of  refraction  n and  the  extinction  coefficient  k,  both  of  which 
may  be  strong  functions  of  wavelength.  These  two  quantities  constitute  the 
real  and  imaginary  parts  of  the  complex  refractive  index 


m 


n + ik 


(11-13) 


To  avoid  confusion  with  extinction  in  the  sense  of  scattering  plus  absorption 
by  small  particles,  the  term  extinction  coefficient  for  k will  not  be  used. 

A negative  sign  is  frequently  found  in  the  complex  refractive  index  relating 
to  the  other  choice  for  the  phase  convention  of  a sinuosidal  wave. 

In  regions  of  low  absorption,  such  as  the  visible  region  for  clear  glasses, 
n and  k are  easily  determined  by  a combination  of  reflectance  and  transmission 
measurements,  or  by  a transmission  measurement  coupled  with  a measurement  of 
minimum  deviation  through  a prism  of  the  material.  In  regions  where  the 
absorption  coefficient  is  very  high,  such  as  near  10  ym  for  the  silicates  of 
concern  in  this  work,  transmission  cannot  be  measured  except  on  thin  films. 

The  most  common  measurement  is  almost-normal  incidence  reflectance  over  a range 
of  wavelengths.  The  data  are  usually  analyzed  to  obtain  the  optical  constants 
by  one  of  the  following  methods: 

(1)  Kramers -Kronig  analysis 

(2)  Multiple  oscillator  fits  to  the  data 

The  Kramers-Kronig  dispersion  equations  relate  the  real  and  imaginary  parts 
of  any  complex  response  function,  under  very  general  conditions  (Stern  1963). 

A response  function  relates  a cause  to  an  effect.  The  comp' ex  amplitude 
reflection  coefficient,  r(w)  = r exp(i  0),  is  such  a suitable  response 
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function  as  is  the  function  1 

*v# 

In  r = In  r(w  ) + i0(w). 

The  Kramers  Kronig  relations  between  the  real  and  imaginary'  parts  of  this 
function  are  given  by 


0(w) 


2“_  p flU 
IT  J bi' 


r(a>‘) 


'2  0)2 


dw" 


(11-14) 


where  P denotes  the  principle  value  of  the  integral.  This  can  be  put  in 
the  following  form  used  in  this  work: 


f] 


CO 

r . 


In  r [X')  -In  r(X)  ^ 

\'2-  x2 


(11-15) 


Since  the  magnitude  of  the  amplitude  reflection  coefficient  is  related  to 
the  measured  intensity  reflection  coefficient  R by 


r = 


(11-16) 


the  6 's  can  be  calculated  if  reflectance  is  known.  In  the  theory  the 
reflectance  must  be  known  from  wavelengths  ofo  to  but  in  practice 
suitable  extrapolations  can  give  a good  approximation  to  the  true  phase 
shifts  0 if  reflectance  measurements  extend  a considerable  distance  on  either 
side  of  the  wavelength  of  interest.  Once  the  phase  shifts  have  been  calculated 
at  a given  X from  the  integral  over  all  wavelengths,  the  optical  constants 
can  be  calculated  at  this  wavelength  from  the  following  expressions: 


n(X) 


1 - r*(  X ) 


1 + r2(x)  - 2r(x)  cos0(x)  (11-17) 


n 


V-'-t  „ . 


The  exact  nature  of  the  extrapolations  used  and  the  computer  program  are 
given  in  Steyer  (1974). 

The  second  method  for  analyzing  reflectance  data  for  optical  constants 
is  to  assume  the  oscillators  responsible  obey  a multiple,  Lorentzian  oscillator 
theory.  This  is  known  to  be  the  case  to  a close  approximation  for  many  simple 
solids  in  the  infrared  lattice  vibration  region.  In  the  Lorentzian  model  the 
complex  dielectric  function  is  given  by 


e*  = 


ic. 


I- 


f 2 

T .0)  • 

J J 


2 2 
-w 
J 


+ itoT- 

J 


(11-19) 


where  the  parameters  (o-,  r^,  and  f.  determine  the  position,  width,  and 

strength  of  the  jth  oscillator,  and  is  the  high  frequency  dielectric  function 
(visible  frequencies  in  this  work).  The  optical  constants  are  related  to  Cj 
and  ^2  by 

e 2 .2 

1 = n ‘ k (11-20) 

e2  = 2nk 


and  the  normal  incidence  reflectance  is  given  by 

R = (n  - l)2  + k2 
(n  + l)2  + k2  . 


(11-21) 


Optical  constants  are  derived  from  the  daia  by  repeatedly  varying  the  oscillator 
parameters  and  calculating  reflectance  spectra  until  a good  fit  is  obtained  to 
the  measured  reflectance  spectrum.  The  oscillator  parameters  so  determined 
produce  the  required  optical  constants. 

In  this  work  the  Kramers -Kronig  technique  was  used,  except  in  the  case  of 
amorphous  olivine  material  produced  by  radiation  damage.  In  this  case  an 
oscillator  fit  was  easier,  since  the  radiation  resulted  in  a thin,  amorphous 
layer  on  top  of  a crystalline  bulk  having  previously-determined  oscillator 
parameters . 
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III.  GLASS  BUBBLES 

The  thin-walled  glass  hubbies  or  micro-balloons  used  in  these 

studies  were  acquired  from  two  different  manufacturers  by  Dr.  Carl 

Dinerman  of  the  Naval  Weapons  Support  Center,  Crane,  Indiana.  They 

are  identified  as  follows: 

E & C:  Emerson  and  Cumming,  Inc. 

Ecospheres  FTF15 

3M  : 3M  Brand 

Glass  Bubbles 

We  will  frequently  use  the  shortened  form  of  the  names  as  given  above. 
The  following  table  summarizes  the  appropriate  information  from  the 
respective  manufacturers'  specification  sheets: 


Table  I 

Properties  of  the  Glass  Bubbles 

E & C 3M 


Reference 

Particle  size  range 

True  Density  (liquid  displacement) 


Tech.  Bull.  14-2-4D 

2 to  44  vim 

.28  to  .32  gm/cm^ 


Technical 
Data  1 


Type  of  material  insoluble  glass  similar  to 

window  glass 

Inner/outer  radii 
(calculated  from  density) 


A.  Size  Distributions 

In  order  to  proceed  with  computer  calculations  of  extinction,  it 
is  necessary  to  have  more  complete  size  information  (distribution  of 
inside  and  outside  radii).  For  this  purpose  both  light  microscopy  and 
scanning  electron  microscopy  were  completed  for  both  samples  of  bubbles. 
Examples  of  the  photomicrographs  are  shown  in  Figs.  1 through  5. 


Figs.  1 and  2 are  light  micrographs  of  the  two  types  of  bubbles.  These 
pictures  have  been  used  to  determine  size  distributions  for  the  particles. 
Measurements  of  individual  particle  diameters  were  fed  into  a simple 
statistics  program  to  determine  the  properties  of  the  size  distributions 
which  are  summarized  in  Fig.  6.  The  mean  sizes  are  essentially  the  same 
for  bubbles  from  the  two  manufacturers  — 19  pm.  The  size  distributions 
summarized  by  the  two  histograms  are  slightly  different.  We  note  a slight 
discrepancy  between  our  mean  size  of  19  pm  and  the  Emerson  and  Cumming 
statement  of  average  size  as  less  than  15  pm,  but  this  difference  is  not 
expected  to  be  important. 

Figs.  4 and  5 show  an  interesting  case  of  a broken  bubble  as  viewed 
by  the  scanning  electron  microscope.  Fig.  5 permits  a rather  accurate 
measurement  to  be  taken  of  the  shell  thickness  in  comparison  to  the  outer 
diameter  crom  this  case  we  have  determined  that  the  ratio  of  inner  to 
outer  radius  is  r./r  = 0.96  which  is  in  excellent  agreement  with  the 
value  determined  from  the  true  density  as  quoted  by  E & C.  Because  of 
the  difficulty  of  determining  shell  thicknesses  and  the  complexity  of 
including  differing  distributions  of  inner  and  outer  radii  we  have  chosen 
to  consider  the  above  ratio  of  r_j/rQ  to  be  representative  of  the  bubbles, 
independent  of  outer  radius  distribution.  This  fact,  along  with  the 
distribution  of  rQ  specifies  the  geometrical  parameters  needed 

for  the  calculations. 

B.  Measurements  of  Optical  Constants 

In  the  vicinity  of  10  pm  all  silicate  materials  have  very  strong 
absorption  bands  due  to  interaction  of  the  light  with  the  stretching 
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Figure  1 

Emerson  and  Cumming  glass  bubbles 
taken  with  the  light  microscope 


Figure  2 

3M  glass  bubbles  taken  with 
the  light  microscope 
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Figure  3 


Emerson  and  Cumming  glass  bubbles  taken  with  the  scanning  electron 

microscope  (^OOx) 


EMERSON  AND  CUMMINGS  GLASS  BUBBLES 


I 

I 


Average  diameter  19.2  g 
Standard  deviation  7.2  p 


3M  GLASS  BUBBLES 


Average  diameter  19.7  g 


Standard  deviation  12. A p 


Figure  6 \ 

Size  distribution  histograms  for  the  two  bubble  materials 

v 

taken  from  photomicrographs. 
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modes  of  vibration  of  the  Si-0  bonds.  This  is  an  example  of  a situation 
where  optical  constants  are  most  easily  derived  from  reflectance  spectra 
measured  at  near-normal  incidence  from  a polished  surface  of  bulk  mate- 
rial and  use  of  the  Kramers -Kronig  analysis.  This  is  the  technique  we 
have  used  to  determine  optical  constants  of  the  two  bubble  glasses. 

Toward  this  end  it  was  necessary  to  melt  the  bubble  material  into 
a homogeneous  mass.  This  apparently  simple  task  proved  to  be  one  of  the 
greatest  challenges  in  the  entire  project.  Simple  heating  of  the  compacted 
material  to  1300  °C  failed  to  cause  them  to  flow  into  anything  like  a 
clear  mass.  After  trying  various  techniques  it  was  found  that  an  oxy- 
hydrogen  torch,  carefu'ly  played  over  a pile  of  the  bubbles  on  a carbon 
block  could  melt  the  material  into  relatively  gas-bubble-free  pieces  of 
about  5 mm  or  more  across.  The  two  bubble  materials  behaved  quite  differ- 
ently under  the  flame,  however.  Gentle  heating  (to  perhaps  600  °C)  of 
the  3M  bubbles  caused  the  material  to  slowly  fuse  and  release  most  of 
the  trapped  gaoses.  The  E & C bubble  material  had  to  be  heated  to  a 
white  hot  heat  (an  estimated  1300  °C)  before  fusion  into  one  molten 
lump  could  be  accomplished.  This  fact  led  us  to  believe  that  the  E & C 
glass  was  probably  quite  similar  to  fused  silica,  although  the  manufactur- 
er only  states  that  it  is  "insoluble  glass". 

In  order  to  provide  further  information  on  the  nature  of  the  glasses, 
refractive  index  (n)  measurements  were  carried  out  at  several  visible 
wavelengths.  The  determination  was  made  by  the  method  of  minimum  devia- 
tion through  a prism  (see  for  example  Jenkins  and  White  1S57).  Two 
flat  faces  were  ground  and  polished  on  samples  of  the  fused  glasses.  The 
angle  of  minimum  deviation  through  each  prism  and  the  prism  angles  were 
measured  with  a mechanical  spectroscope  that  could  be  read  to  ^ minute  of 
an  arc  (.008°),  fitted  with  a crystal -orientating  goniometer.  Results 
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calculated  from  these  data  for  n at  several  wavelengths  in  th''  visible 
are  shown  in  Table  II  below. 

Table  II 

Measured  Refractive  Indices  in  the  Visible  for  the  Two 
Bubble  Materials  and  for  Fused  Silica 


Wavelength 

Fused  Silica 

E & C 

3M 

.63  mp 

1.4572 

1.4572  * .0005 

1.504 

+ .005 

.59 

1.4584 

1.4586 

.56 

1 .4596 

1.4596 

1.494 

+ .01 

.44 

1.4664 

1.4660 

1.508 

i .01 

The  larger  experimental  uncertainties  in  the  case  of  the  3M  material  was 
due  to  the  lack  of  transparency  in  this  sample  because  of  included  gas 
bubbles.  It  is  apparent  from  the  results  that  the  E & C material  after 
being  fused  by  our  technique  is  optically  similar  to  fused  silica  within 
narrow  limits.  This  confirms  the  chemical  similarity  to  silica  already 
observed  in  the  melting  studies.  On  this  basis  the  early  progress  reports 
of  this  work  used  published  optical  constants  in  the  10  pm  region  for  silica  in 
calculating  extinction  of  bubble  material.  The  3M  material  after  being 
fused  by  our  technique  is  quite  different  in  refractive  index,  being 
fairly  near  to  values  quoted  for  borosii !cate  crown  glass. 

Reflectance  m°asurements  on  polished  faces  of  the  two  fused  bubble 
materials  were  measured  from  3 pm  to  25  pm.  The  almost-normal  incidence 
measurements  were  taken  with  two  specular  reflectance  attachments  fitted 
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into  a Perkin  Elmer  337  and  a Beckman  IR12  spectrophotometer.  An  alumi- 
num reflectance  standard  was  used  in  a sample-in-standard-in  method. 
Resulting  reflectance  spectra  in  the  7pn  to  25  pm  region  for  the  two 
materials  and  for  fused  silica  are  presented  in  Fig.  7.  Again  one  ob- 
serves that  the  E & C material  is  quite  similar  to  fused  silica  except 
for  the  magnitudes  of  the  reflectance  peaks,  while  the  3M  material  is 
considerably  different.  Complex  indices  of  refraction  derived  from  the 
reflectance  measurements  by  the  Kramers -Kronig  computer  program  are 
summarized  in  Table  Iil. 


Table  III 

Complex  Optical  Constants  for  the  Two  Bubble  Materials 
and  for  Fused  Silica 


Fused 

Silica 

E & C Material 

3M  Material 

X 

£ 

k 

n 

k 

n 

k 

7.8 

.64 

.09 

.88 

.18 

8.0 

.48 

.18 

.46 

.39 

.75 

.25 

8.2 

.38 

.55 

.48 

.55 

.68 

.41 

8.4 

.46 

.80 

.49 

.73 

.64 

.59 

8.6 

.45 

.94 

.49 

.95 

.65 

.77 

8.8 

.39 

1.37 

.50 

1.16 

.65 

.93 

9.0 

.68 

2.31 

.69 

1.63 

.75 

1.25 

9.2 

1.72 

2.51 

1.07 

1.82 

..94 

1.41 

9.4 

2.73 

1.71 

1.51 

1.83 

1.31 

1.59 

9.6 

2.89 

1.03 

2.09 

1.55 

1.61 

1.59 

9.8 

2.76 

.54 

2.25 

1.20 

1.84 

1.59 

10.0 

2.40 

.22 

2.24 

.93 

2.13 

1.24 

10.4 

2.01 

.14 

2.15 

.57 

2.21 

.82 

10.8 

1.82 

.10 

1.99 

.43 

2.17 

.61 

11.2 

1.72 

.14 

1.89 

.35 

11.4 

1.69 

.16 

1.80 

.33 

2.09 

.39 

12.1 

1.66 

.26 

1.66 

.37 

1.84 

.30 

These  constitute  the  optical  parameters  required  as  input  for  the 
extinction  calculations. 


C.  Extinction  Calculations 


Using  the  two  sets  of  measured  bubble  parameters  (size  distri- 
butions and  complex  refractive  indices)  presented  in  the  previous  two 
sections  as  input  for  the  concentric  spheres  computer  program,  extinction 
was  calculated.  Size  distributions  used  were  those  of  Fig.  6 for  outer 
radii,  with  the  inner  radii  set  by  the  assumption  that  r./r  = 0.96, 
and  optical  constants  us.  d were  from  Table  III.  Results  are  shown' in 
Fig.  8.  The  particular  wavelengths  of  interest  are  those  in  the  9.6  and 
10.6  ym  COg  laser  bands  sketched  in  at  the  bottom  of  the  curve  from  one 
laser  manufacturer's  specifications.  Generally  speaking  the  volume 
normalized  extinction  cross  sections  are  0.8  to  1.0  x 10^  cm~^  in  both 
bands,  for  both  materials.  In  some  sense  Fig.  8 summarizes  the  final 
conclusions  of  this  work  as  it  regards  the  specific  glass  bubble  samples 
at  our  disposal.  However,  some  further  insight  can  be  obtained  into  the 
mechanisms  responsible  for  the  extinction  by  avoiding  the  complexity  of 
the  complete  size  distribution  to  examine  extinction  by  different  single 
sizes.  This  insight  may  provide  the  basis  for  optimizing  the  extinction 
efficiency. 

An  important  question  in  this  connection  is  how  does  the  extinction 
efficiency  vary  with  bubble  size.  To  illustrate  this,  we  have  made  a 
series  of  calculations  at  the  fixed  wavelength  of  10  ym  (near  the  maximum 
extinction  in  both  materials)  for  a series  of  increasing  bubble  radii,  the 
relative  shell  thickness  being  kept  constant  as  previously  (r./r  = 0.96). 
The  results  are  shown  in  Fig.  9.  One  observes,  perhaps  surprisingly,  that 
the  smaller  bubbles  are  the  more  effective.  There  is  no  peak  in  the 


Figure  9 

Extinction  normalized  per  unit  volume  of  solid  calculated  for 
different  outside  radii  of  3M  bubbles  with 
r{/r„  « 0.96. 
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extinction  efficiency  as  the  size  approaches  the  wavelength  as  there  would 
be  if  the  optical  constants  were  indeed  constant,  and  if  k were  small. 

The  fact  is  that  the  large  variations  in  optical  constants  in  the  Si-0 
resonance  region  produce  strong  absorption  which  is  more  important  than 
the  scattering.  Table  IV  shows  how  the  extinction  efficiency  factor  is 
distributed  between  scattering  and  absorption  efficiencies. 


Table  IV 


Scattering  and  Absorption 

Contributions  to 

Extinction 

for  3M  Glass  Bubbles  (r  = 5 ym,  r 

•-  = 4.8  ym) 

X(pm) 

^ext 

^sca 

^abs 

9.1 

.76 

.17 

.59 

9.3 

.96 

.22 

.73 

9.6 

1.21 

.30 

.91 

10.0 

1.29 

.35 

.94 

10.5 

1.08 

.33 

.75 

11.1 

.77 

.25 

.51 

One  sees  that  absorption  is  the  more  important  contribution  to  extinction 
even  though  the  particle  size  is  comparable  to  the  wavelength.  Furthermore, 
since  absorption  is  more  important,  the  extinction  per  unit  volume  of 
glass  increases  as  the  effective  thickness  of  the  glass  decreases.  This 
is  an  example  of  a type  of  absorption  saturation  in  this  strongly  absorb- 
ing region  (see  Day,  et  al,  1974  for  a discussion  of  this  saturation). 

Once  the  particle  absorbs  essentially  all  the  radiation  falling  on  it 
(speaking  rather  loosely  for  simplicity)  further  increases  in  thickness 
will  cause  little  increase  in  extinction,  though  the  volume  (and  mass) 
will  increase  substantially.  Thus  the  volume-  or  mass-  normalized  extinc- 
tion must  decrease  as  the  effective- absorbing  thickness  increases.  The 


19 


conclusion  must  follow  that  the  smaller  bubbles,  having  smaller  wall 
thicknesses,  are  the  more  effective.  An  obvious  improvement  in  extinction 
could  therefore  be  affected  by  using  smaller  size  fractions  of  the 
bubbles.  With  this  in  mind  we  have  calculated  the  extinction  spectrum 
for  5 pm  radius  bubbles  (again  with  r../ro  = 0.96)  with  optical  constants 
appropriate  to  fused  silica,  3M  material,  and  E & C material.  Fig.  10 
presents  the  results.  One  notices  the  expected  increase  of  maximum 
by  a factor  of  about  2.  Of  course  this  increased  extinction  efficiency 
for  silicate  glass  materials  theoretically  could  be  enhanced  even  more 
by  dispersing  the  glass  as  very  fine  particles.  To  illustrate  this  we 
have  included  in  Fig.  10  the  results  of  one  Hie  calculation  (which  is 
practically  identical  to  a Rayleigh  calculation  for  this  size)  for  .02  pm 
radius  E & C spheres.  There  is  a further  enhancement  of  extinction,  al- 
though the  peak  has  shifted  to  about  8.8  pm.  This  case  illustrates  the 
maximum,  volume  normalized  extinction  that  could  be  expected  for  any  kind 
of  particles  made  from  the  3M  bubble  material.  In  such  a hypothetical 
case,  however,  we  have  gotten  away  from  the  original  idea  of  scattering 
glass  bubbles  and  are  considering  a smoke  of  silicate  glass  such  as  will 
be  considered  in  the  next  section. 
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Figure  10 

Extinction  normalized  per  unit  volume  of  the  solid  calculated 
for  5 pm  bubbles  of  the  three  materials.  Also 
shown  is  a calculation  for  .02pm  solid  spheres 
of  the  E & C material. 
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IV.  SILICATE  SMOKE 

In  the  proposal  for  this  contract  the  author  suggested  very 

small  particles  (.01  to  0.5  ym)  of  silicate  materials,  called  here 

silicate  smokes,  as  effective  agents  for  extinction  of  CC^  laser 

radiation.  This  suggestion  was  based  on  previous  work  done  in  this 

laboratory  in  investigation  of  possible  interstellar  grain  materials 

which  have  a strong  extinction  "bump"  in  the  10  ym  spectral  region. 

This  has  been  attributed  to  absorption  caused  by  the  Si-0  vibration 

modes  in  silicates.  (For  a review  of  interstellar  grains  and  optical 

studies  of  small  particles  see  Huffman  1977).  Theoretical  extinction 

based  on  measured  optical  constants  had  shown  volume  extinction  coef- 

5 -1 

ficients  approaching  10  cm  . The  particular  one  of  many  possible 
silicates  suggested  was  the  mineral  olivine,  (Mg, Fe) ^Si 0^ » for  the 
following  two  reasons: 

(1)  Olivine  re-condenses  without  decomposition  after  being 
vaporized  to  form  smoke. 

(2)  It  has  several  major  absorption  peaks  near  9.6  and  10.6ym. 

After  the  work  got  under  way  it  became  apparent  that  the  olivine 

smoke  sometimes  condenses  as  crystalline  particles  and  other  times  as 
amorphous  particles.  The  two  forms  have  significantly  different  infra- 
red properties.  Although  the  optical  constants  of  crystalline  olivine 
had  already  been  determined  in  this  laboratory  (Steyer  1974),  amorphous 
olivine  had  not  been  measured,  and  there  were  serious  difficulties  stand- 
ing in  the  way  of  doing  it.  Apparently  no  one  had  been  able  to  produce 
large  (cm-sized)  chunks  of  amorphous  olivine.  As  do  most  solids,  olivine 
much  prefers  to  solidify  in  a crystalline  form.  The  solution  to  the 
sample  preparation  problem  came  through  another  project  in  which  we  were 
exposing  olivine  crystals  to  high  radiation  doses  of  heavy  ions  from 


a Van  de  Graaff  accelerator,  in  an  effort  to  produce  astrophysically 
related  radiation  damage  effects.  The  high  degree  of  radiation  damage 
produced  an  essentially  amorphous  layer  on  the  bulk  crystal.  Measure- 
ments of  optical  constants  were  therefore  conducted  on  the  amorphous 
olivine  to  enable  extinction  calculations  to  be  completed  on  this  kind 
of  particle. 

The  olivine  smoke  particles  were  produced  by  vaporizing  mm-sized 
pieces  of  natural  olivine  from  the  Globe,  Arizona,  area.  The  crystal 
fragments  are  embedded  in  hollowed-out  carbon  arc  electrodes,  and  an  ac 
or  dc  arc  is  struck  in  air  at  a current  of  about  10  amp.  The  smoke 
produced  as  the  vaporized  olivine  condenses  was  usually  caught  on  infra- 
red-transparent  substrates  of  KC1  or  silicon  for  transmission  measure- 
ments and  for  mass  determinations  using  an  analytical  balance.  From 
these  results  extinction  ( c^)  could  be  determined.  Samples  for  electron 
microscopy  were  occasionally  taken  by  placing  standard  electron  micro- 
scope grids  in  the  smoke  stream.  X-ray  diffraction  patterns  were  taken 
on  the  collected  material  using  a powder  camera.  This  revealed  that  the 
particles  were  usually  highly  crystalline,  with  identical  lattice  spacings 
to  the  crystal  feed  material.  Occasionally,  however,  particles  were  pro- 
duced that  showed  no  x-ray  diffraction  pattern  at  all,  and  were  thus  con- 
cluded to  be  quite  amorphous.  The  conditions  for  production  of  this 
amorphous  material  were  never  clearly  determined,  but  appeared  to  be  re- 
lated to  the  temperatures  in  the  condensation  zone.  Factors  such  as  the 
shape  of  the  electrode  tips  and  the  exact  position  of  the  molten  ball  of 
material  appeared  to  be  important  variables,  but  these  were  very  diffi- 
cult to  control  in  the  somewhat  violent  process.  Infrared  transmission 
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spectra  of  the  two  smoke  forms  were  always  distinctively  different  in 
the  10  vm  region.  The  crystalline  smoke  showed  at  least  two  well  defined 
absorption  peaks  while  the  amorphous  smoke  showed  only  one  broad  peak. 

The  infrared  transmission  measurements  taken  on  the  smoke-coated 
substrates,  plus  mass  determinations  made  by  weighing  the  coated  and 
uncoated  substrates,  permitted  extinction  coefficients  to  be  calculated. 
Extinction  calculations  and  measurements  were  accomplished  for  both 
crystalline  and  amorphous  olivine  smokes.  These  results  are  presented 
in  the  following  section. 

A.  Size  Distributions 

Size  distributions  for  the  olivine  smoke,  necessary  for  extinction 
calculations,  had  been  made  in  earlier  studies  at  this  laboratory  (Day 
1975)  from  electron  micrographs.  The  particles  are  spherical,  ranging  in 
size  from  about  .01  pm  to  about  0.2  pm  radius  with  mean  radius  of  about 
.03  pm.  Since  the  volume  normalized  extinction  coefficients  are  inde- 
pendent of  size  in  the  small  particle,  Rayleigh  limit,  the  exact  size 
distribution  used  in  the  calculations  is  immaterial  and  is  not  given  in 
this  report. 

B.  Optical  Constants  Determinations 

Because  of  the  crystalline  anisotropy  of  olivine  it  is  necessary 
that  three  sets  of  optical  constants  be  determined  — for  the  electric 
vector  of  the  light  parallel  to  each  of  the  three  crystal  axes.  This 
had  been  done  previously  at  this  laboratory  (Steyer  1974)  by  Kramers- 
Kronig  analysing  spectral  reflectance  data  taken  with  an  x-ray  oriented 
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crystal  oriented  in  the  three  different  ways  with  respect  to  a polarizer 
in  the  reflectometer.  Results  for  the  optical  constants  of  the  crystalline 
material  are  shown  in  Fig.  11  in  terms  of  the  complex  dielectric  functions 
e-j  and  ^ The  dielectric  functions  are  used  here  because  of  the  simpli- 
city of  the  Rayleigh  extinction  equation  written  in  terms  of  them. 

In  the  case  of  the  amorphous  material,  optical  constants  were 
determined  from  the  reflectance  spectra  of  the  highly  disordered  surface 


layer  produced  on  a polished  piece  of  single  crystal  olivine  by  ion  bom- 

06 


bardment.  The  polished  crystal  .was  subjected  to  a dose  of  about  5 x 10 


2 + 

ions/cm  of  Ne  accelerated  to  1.5  Mev  in  a Van  de  Graaff  accelerator 


of  the  Physics  Department,  University  of  Arizona.  Interference  fringes 
in  visible  light  determined  that  the  damaged  layer  was  about  2 ym  thick, 
in  good  agreement  with  values  calculated  from  the  expected  range  of 
such  atomic  particles.  With  our  prior  knowledge  of  the  oscillator  param- 
eters necessary  to  describe  the  undamaged  underlying  layer,  and  the 
thickness,  oscillator  parameters  for  the  damaged  layer  were  varied  until 
a best  fit  was  achieved  for  the  reflecting  surface.  Fig.  12  shows  the 
measured  reflectance  of  the  undamaged  and  the  damaged  surface,  along  with 
a calculated  reflectance  for  an  infinitely  thick  layer  of  damaged  mater- 
ial. These  data  are  for  one  polarization  direction  of  the  light  inci- 
dent on  the  sample.  Resulting  optical  dielectric  functions  are  shown 
in  Fig.  13  resulting  from  the  best  fit  oscillator  parameters.  Measure- 
ments similar  to  those  in  Fig.  12  (middle)  were  made  for  the  other  two 
orientations  of  the  major  axes  and  optical  constants  derived  from  the 
best  fits.  The  error  bars  in  Fig.  13  represent  the  spread  in  values 


of  e-j  and  obtained  for  the  three  different  determinations.  This  work 


is  being  submitted  for  publication  (Kratschmer  and  Huffman  1978} 


DIELECTRIC  FUNCTIONS 
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The  dielectric  functions  for  single  crystal  olivine  of  Fig.  11  and  for 
amorphous  olivine  of  Fig.  13  are  those  used  in  calculations  of  extinction 
presented  in  the  next  section. 


C.  Extinction  Measurements  and  Calculations 

With  the  optical  constants  for  crystalline  and  amorphous  olivine 
taken  from  the  preceeding  section,  extinction  calculations  have  been  made 
for  the  two  kinds  of  smoke.  The  results  are  shown  as  the  solid  lines  in 
Figs.  14  and  15.  The  results  of  measurements  of  volume  normalized  extinc- 
tion coefficients  for  the  two  smokes  also  are  displayed  in  Figs.  14  and 
15  by  the  dotted  lines.  These  two  figures  summarize  the  main  results  of 
this  study  as  it  relates  to.  the  silicate  smokes  Positions  of  the  9.6 
and  10.6  ym  CO^  laser  bands  are  shown  on  the  fit  ~es  as  previously.  The 
appearance  of  the  extinction  curves  for  crystalline  and  amorphous  mater- 
ials are  seen  to  be  distinctly  different.  There  are  three  major  peaks 
in  the  calculated  extinction  for  crystalline  smoke  and  one  broad  band 
for  the  amorphous  smoke. 

Some  helpful  insights  into  the  nature  of  these  extinction  bands  can 
be  obtained  by  looking  at  their  cause  on  the  basis  of  the  extinction  equa- 
tion for  ellipsoidal  particles  in  the  Rayleigh  limit  (r«X,  nr«A  ). 

The  nice  thing  about  the  Rayleigh  equation  is  that  it  is  simple  enough  that 
one  can  see  the  effect  of  varying  its  parameters,  and  it  is  not  restricted 
to  spheres,  so  that  insights  into  shape  effects  are  possible.  The  extinc- 
tion cross  section  per  unit  volume  for  an  ellipsoidal  particle  can  be 
written  in  terms  of  the  complex  dielectric  function  of  the  solid  and  the 
quantities  4ttL.,  known  as  depolarization  factors  in  electromagnetic 
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theory. 
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{Lj£l  + (l-L^.)}2  + (Lj£2)2 


(IV-1) 


The  above  assumes  that  the  particle  is  in  free  space  and  that  scattering 

is  negligible.  For  a general  ellipsoid  there  are  three  different  L.'s 

J 

which  are  subject  to  the  conditions  0<L . <1  and  Z L.  =1  {see  van 

J J 

de  Hulst  1957),  related  to  the  lengths  of  the  semi-axes.  One  sees  that 
a resonance  occurs  in  the  extinction  when 


Lj.e,  + (1  - L,)  - 0, 


(IV-2) 


the  strength  of  the  resonance  being  stronger  as  is  smaller.  For 

spheres  for  which  L.  = 1/3  for  all  the  axes,  the  resonance  condition  is 

J 


e,  - -2 


These  resonances,  which  occur  at  the  wavelengths  for  which  e-j  = -2 
for  spheres  and  which  occur  in  the  negative  e-j  region  for  all  shapes, 
have  sometimes  been  called  surface  modes.  They  are  probably  more  correct- 
ly termed  Fr'Ohl ich  modes  (after  Frohlich  1937),  A review  with  examples 
is  given  in  Huffman  (1977).  These  modes  occur  in  the  vicinity  of,  but 
not  exactly  at,  the  position  of  the  strong  infrared  lattice  mode  reso- 
nances. The  condition  that  the  strength  goes  up  as  c2  at  the  Frohlich 
mode  frequency  goes  down  means  that  extremely  high  extinction  resonances 
could  occur. 

Returning  now  to  the  smoke  extinction  curves  of  Fig.  14,  the  three 
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peaks  in  the  extinction  calculations  occur  at  precisely  the  wavelengths 
where  = -2  for  the  three  polarization  directions  of  the  crystal. 

This  could  be  determined  by  a careful  comparison  of  Figs.  11  and  14, 
but  is  more  clearly  shown  in  Steyer's  thesis  (1974).  The  calculated 
extinction  peak  of  the  amorphous  smoke  (Fig.  15)  falls  at  the  wave- 
length where  is  minimum,  although  it  never  quite  reaches  -2. 

Because  the  resonances  in  z^  are  narrower  in  the  crystalline  than  they 
are  in  the  amorphous  olivine,  the  e^'s  at  wavelengths  where  e-j  = -2 
are  less,  giving  higher  peak  efficiencies.  Going  along  with  the  high- 
er peak  efficiencies,  however,  is  the  fact  that  the  extinction  resonances 
are  much  narrower.  The  appearance  of  the  two  extinction  spectra  are  thus 
considerably  different. 

In  both  cases  the  peak  extinction  efficiencies  are  less  by  a factor 
of  2 to  4 in  the  measurements  as  compared  to  the  calculations.  This  seems 
to  have  been  a general  trend  in  the  few,  previously  reported  cases  in 
which  attempts  have  been  made  to  compare  theoretical  and  experimental 
extinctions  with  no  arbitrary  scaling  (see  Steyer  et  al  1974,  and  Day  1975). 
The  explanation  for  the  discrepancy  seems  to  lie  in  the  breakdown  of  the 
independent  sphere  approximation  for  smoke  loosely  trapped  on  substrates. 
Again,  Eqs.  IV-1  and  2 give  some  insight  into  what  might  happen  if  two 
or  more  particles  stick  together.  They  would  act  somewhat  like  an  elon- 
gated particle  approximated  by  different  L-'s,  giving  Frohlich  modes 

J 

displaced  from  the  e-j  = -2  wavelength,  but  still  in  the  negative 
e-|  region.  A distribution  of  shapes,  even  if  caused  by  clumping  of  spheres, 
should  give  rise  to  a broadened  extinction  band  whose  maximum  is  now 
depressed.  This  is  exactly  what  one  sees  in  comparing  the  sphere 
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calculations  and  the  measurements  of  Figs.  14  and  15.  Even  before  the 
smoke  particles  are  trapped  on  the  substrate  they  form  themselves  into 
chains  like  strings  of  pearls.  This  appears  to  happen  in  the  near  vi- 
cinity of  the  arc  and  prevents  the  spherical  particles  from  acting 
independently.  An  aerosol  cf  the  smoke  might  have  slightly  higher  ex- 
tinction than  when  coating  a substrate;  however,  a good  guess  is  that 
an  aerosol  of  the  olivine  smoke  particles  would  show  extinctions  somewhere 
between  the  measured  and  calculated  curves,  when  first  produced.  After 
long  periods  of  time, agglomeration  of  the  smoke  would  likely  decrease 
the  volume  extinction  coefficient  even  more. 


The  highest  extinctions  are  seen  in  the  calculations  for  crystalline 


smoke,  which  reach  to  4.4  x 10^  cm*^  at  11.15  pm.  Although  the  peaks 


of  the  two  main  absorptions  are  near  to  9.6  and  10.6  pm  bands,  they 
are  not  exactly  centered  on  the  laser  bands.  Due  to  the  sharpness  of  the 
absorption  bands  their  slight  displacement  from  the  desired  values  pro- 
duces CO^  laser  extinction  values  considerably  below  their  peak  values. 
The  breadth  of  the  extinction  band  for  amorphous  smoke  causes  it  to  be 
almost  equally  effective  in  both  of  the  CO,,  bands.  Both  sets  cf  measure- 
ments are  accidentally  about  the  same  at  9.6  and  10.6  pm. 


V.  COMPARISONS  AND  CONCLUSIONS 


■ • Having  presented  in  the  above  sections  the  results  of  extinction 
for  the  two  bubble  materials  and  the  two  silicate  smokes,  it  is  now 
possible  to  give  a more  educated  assessment  cf  which  material  would  pro- 
vide the  more  effective  extinction  agent.  The  appropriate  results  were 
given  in  Fig.  8 for  the  bubble  materials  and  Figs.  14  and  15  for  the 
silicate  smokes.  Amazingly  and  disgustingly,  all  four  candidate^  give 
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about  the  same  extinction  within  a factor  of  two.  Thus  none  of  the 
materials  are  clear  choices  over  the  others.  Any  choice  among  these 
would  need  to  be  made  on  other  grounds,  such  as  which  material  is  most 
easy  to  produce  and  disperse  in  practical  situations  and  which  aerosol 
resists  agglomeration  and  settling  best. 

A significant  point  regarding  the  bubbles  is  that  scattering 
is  not  the  dominant  extinction  mechanism,  so  that  smaller  sized  fractions 
would  give  more  extinction  per  unit  volume  of  solid.  It  should  also  be 
kept  in  mind  that  the  extinction  per  unit  volume-of-bubbles  would  be 
less  than  the  extinction  per  unit  volume  of  solid  (c*p)  discussed  in  this 
work.  Hence,  if  a pyrotechnic  device  had  to  carry  the  bubbles  already 
made,  as  contrasted  with  materials  to  make  the  bubbles. from,  the  appro- 
priate volume  extinction  coefficient  would  be  less.  If  on  the  other  hand 
the  weight  is  the  most  important  parameter  in  the  aerosol  producer,  it 
would  not  matter  whether  the  device  carried  pre-made  bubbles  or  produced 
bubbles  from  bulk  material. 

In  view  of  the  similarity  of  extinction  coefficients  derived  it 
appears  that  other  considerations  must  distinguish  between  the  use  of  one 
or  the  other  of  these  materials.  The  absorption  mechanism  responsible  for 
the  strong  extinction  resonances  in  the  silicate  smokes  could,  however,  be 
higher  in  other  materials.  Complex  optical  constants  are  required 
before  the  determination  of  effectiveness  for  a given  solid  can  be  made. 
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